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Abstract 

We report on detailed Monte Carlo comparisons of selection variables to separate tbH^ signal events from 
the Standard Model ti background using an Iterative Discriminant Analysis (IDA) method. While kinematic 
differences exist between the two processes whenever m H ± ^ m w ± , the exploration of the spin difference 
between the charged Higgs and the gauge boson becomes crucial in the particularly challenging case of 
near degeneracy of the charged Higgs boson mass with the mass. The TAUOLA package is used to decay 
the tau leptons emerging from the charged Higgs and boson decays taking the spin difference properly into 
account. We demonstrate that, even if the individual selection variables have limited discriminant power, the 
IDA method achieves a significant separation between the expected signal and background. For both Tevatron 
and LHC energies, the impact of the spin effects and H ± mass on the separation of signal and background has 
been studied quantitatively. The effect of a hard transverse momentum cut to remove QCD background has 
been studied and it is found that the spin effects remain important. The separation is expressed in purity versus 
efficiency curves. The study is performed for charged Higgs boson masses between the W mass and near the 
top mass. 
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1. INTRODUCTION 

The importance of charged Higgs boson searches has in the recent years been emphasized Hl|2l[3jffl for 
LEP, a future International Linear Collider (ILC), the Tevatron and the Large Hadron Collider (LHC), 
as the detection of a charged Higgs boson would be a definite signal for the existence of New Physics 
beyond the Standard Model (SM). Charged Higgs bosons naturally arise in non-minimal Higgs scenarios, 
such as Two-Higgs Doublet Models (2HDMs). A Supersymmetric version of the latter is the Minimal 
Supersymmetric Standard Model (MSSM). It is a Type II 2HDM with specific relations among neutral 
and charged Higgs boson masses and couplings, dictated by Supersymmetry (SUSY) Q. 

The Tevatron collider at Fermilab is currently in its second stage of operation, so-called Run 2, 
with a center-of-mass (CM) energy of y/s = 1.96 TeV. This machine will be the first one to directly 
probe charged Higgs boson masses in the mass range up to m H ± ~ m t . Starting from 2008, the LHC at 
CERN will be in a position to confirm or rule out the existence of such a particle over a very large portion 
of both the 2HDM and MSSM parameter space, m H ± < 400 GeV, depending on tan j3, the ratio of the 
vacuum expectation values of the two Higgs doublets (see the reviews (6l|71[H anc ^ a recent study 0). 

At present, a lower bound on the charged Higgs boson mass exists from LEP ifTOl . m#± > m\y±, 
independently of the charged Higgs boson decay Branching Ratios (BRs). This limit is valid within any 
Type II 2HDM whereas, in the low tan (3 region (below about 3), an indirect lower limit on m#± can 
be derived in the MSSM from the one on itla (the mass of the pseudoscalar Higgs state of the model): 
m 2 H± « m^ v± + m\ > (130 GeV) 2 . 

If the charged Higgs boson mass m H ± satisfies m H ± < m t — mi,, where m< is the top quark mass 
and m& the bottom quark mass, bosons could be produced in the decay of on-shell (i.e., T t — > 0) top 
(anti-)quarks t — > bH + , the latter being in turn produced in pairs via gg fusion and qq annihilation. This 
approximation is the one customarily used in event generators when m H ± < m t . Throughout this study 
we adopt the same notation as in Ref . ifTTTl : charged Higgs production is denoted by qq, gg — > tt — > tbH^ 
if due to (anti-)top decays and by qq, gg — > tbH^ if further production diagrams are included. In 
fact, owing to the large top decay width (Tt ^ 1.5 GeV) and due to the additional diagrams which do 
not proceed via direct tt production |[T2l[T3l[T4l . charged Higgs bosons could also be produced at and 
beyond the kinematic top decay threshold. The importance of these effects in the so-called 'threshold' 
or 'transition' region (m H ± « m t ) was emphasized in Les Houches proceedings lfT5l[T6ll as well as in 
Refs. | fTTl[T7l[T8l[T9ll . so that the calculations of Refs. | fi"2l[L3l (based on the appropriate qq, gg -> tbH* 1 
description) are now implemented in HERWIG H3 HU US HI and PYTHIA US HH. A comparison 
between the two generators was carried out in Ref. ifTTIl . For any realistic simulation of production 
with m H ± > m t the use of these implementations is important. In addition, in the mass region near the 
top quark mass, a matching of the calculations for the qq, gg — > tbH^ and gb — > tH^ processes might 
be required ll25ll . 

A charged Higgs boson with m H ± < m t decays predominantly into a r lepton and a neutrino. 
For large values of tan j3 ( > 5) the corresponding BR is near 100%. For m#± > mt, — > TV T IS 
overtaken by —>■ tb, but the latter is much harder to disentangle from background than the former. 
The associated top quark decays predominantly into a boson, or at times a second charged Higgs 
boson, and a b quark. The reaction 

qq^gg^tbH* (t - bW) (f^ - r^) (1) 

is then a promising channel to search for a charged Higgs boson at both the Tevatron (where the dominant 
production mode is qq) and the LHC (where gg is the leading subprocess). If the — * tv t decay 
channel is used to search for Higgs bosons, then a key ingredient in the signal selection process should 
be the exploitation of decay distributions that are sensitive to the spin nature of the particle yielding the 
r lepton (H^ in the signal or in the background), as advocated in Refs. ll26l l27l l28l l29ll (see also 
QUI |3"TT0 . The r spin information affects both the energy and the angular distribution of the r decay 
products. 

In the search for a charged Higgs boson signal containing a r lepton, not only the magnitude of the 
production cross section is important, but also the efficiency of identifying the r lepton in the hadronic 
environment plays a crucial role. Since r leptons have a very short life-time (~ 10~ 6 s), they decay 
within the detectors and can only be identified through their decay products. In about 35% of the cases 
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they decay leptonically and about 65% of the times they do so hadronically. Both of these decay modes 
are usually addressed in charged Higgs boson searches by employing dedicated r lepton triggers. The 
identification of taus in hadronic pp collisions has recently been studied, e.g. Z — > t + t~ events 11321 
and further details are given in fl33"Tl . 

It is the purpose of this note to outline the possible improvements that can be achieved at the Teva- 
tron and LHC in the search for charged Higgs bosons, with mass below the top mass and including the 
appropriate description of the spin effects in the — ► rv T decay. In order to quantify the spin effect an 
Iterative Discriminant Analysis (IDA) method has been applied, which is a powerful tool to separate sig- 
nal and background, even in cases such as the one presently under study when several selection variables 
with limited discriminant power are present. 

2. TEVATRON ENERGY 

We start by studying charged Higgs production qq, gg — > tbH ± with subsequent decays t — > bW, — ► 
tv t at the FNAL Tevatron with yfs = 1.96 TeV. In the following we analyze hadronic decays of the 
boson and r lepton (W^ — ► qq', r — > hadrons+i^ r ), which results in the signature 26+2j+Tj e t+p™ ss (2 
b jets, 2 light jets, 1 r jet and missing transverse momentum). The most important irreducible background 
process is qq,gg — > tt with the subsequent decays t — ► bW + and t — ► b~W~, one boson decaying 
hadronically (W^ — ► gg') and one leptonically (W T — ► rz^ T ), which results in the same final state 
particles as for the expected signal. 

2.1 Simulation and Detector Response 

The signal process qq, gg — ► t6f^ ± is simulated with PYTHIA ll24l . The subsequent decays i — > My 1 * 1 
(or its charge conjugate), — > gff and -fr^ — > r^ r are also carried out within PYTHIA, whereas 
the r leptons are decayed externally with the program TAUOLA Il34ll35l . which includes the complete 
spin structure of the r decay. The background process qq, gg — » tt is also simulated with PYTHIA with 
the built-in subroutines for ti production. The decays of the top quarks and bosons are performed 
within PYTHIA and that of the r lepton within TAUOLA. 

The momenta of the final b and light quarks from the PYTHIA event record are taken as the mo- 
menta of the corresponding jet, whereas for the r jet the sum of all non-lep tonic final state particles as 
given by TAUOLA is used. The energy resolution of the detector and parton shower and hadronization 
effects are emulated through a Gaussian smearing (A(p t ) / Pt) 2 = (O.SO/y 7 ^) 2 of the transverse mo- 
mentum p t for all jets in the final state, including the r jet 0. As typical for fast simulation studies, 
no effects of underlying events, are simulated. Events are removed which contain jets with less than 20 
GeV transverse momentuirH, corresponding to about \rj\ > 3. The transverse momentum of the leading 
charged pion in the r jet is assumed to be measured in the tracker independently of the transverse mo- 
mentum of the r jet. The identification and momentum measurement of the pion is important to fully 
exploit the r spin information. In order to take into account the tracker performance we apply Gaussian 
smearing on l/pl with a{l/pJ)[TeV- 1 ) = v/0.52 2 + 22 2 /(pJ r [GeVJ) 2 sin 9 n , where V is the polar 
angle of the tt. The missing transverse momentum p™ lss is constructed from the transverse momenta of 
all visible jets (including the visible r decay products) after taking the modelling of the detector into 
account. The generic detector description is a good approximation for both Tevatron experiments, CDF 
and DO. 

2.2 Expected Rates 

For completeness we present a brief discussion of the expected cross section of the charged Higgs boson 
signature under investigation. The signal cross section has been calculated for tan (3 = 30 and m H ± = 
80, 100, 130 and 150 GeV with PYTHIA, version 6.325, using the implementation described in ll25l . in 
order to take the effects in the transition region into account. Furthermore, it has been shown in ifTTI that 
the signal cross section for tbH^ agrees with the one from the top-decay approximation tt — > tbH^ for 
charged Higgs boson masses up to about 160 GeV if the same factorization and renormalization scales 
are used. Thus, we have used everywhere in this study the factorization scale (mt + ni}j± ) / 4 and the 

2 In order to be largely independent of the specific detector performance, no requirement on the jet resolution is applied. 



2 



renormalization scale m H ± for both signal and background (i.e., those recommended in |[25l as most 
appropriate for the tbH^ signalJH, since the primary purpose of our study is to single out variables that 
show a difference between our and data samples and that this can unambiguously be ascribed to 
the different nature of the two kinds of bosons (chiefly, their different mass and spin state). In addition, 
the running b quark mass entering in the Yukawa coupling of the signal has been evaluated at m H ± . This 
procedure eventually results in a dependence of our background calculations on tan (5 and, especially, 
m H ± that is more marked than the one that would more naturally arise as only due to indirect effects 
through the top decay width. Hence, the cross sections have been rescaled with a common factor such 
that the total ti cross section is <7^ rod = 5.2 pb ll36l . To be more specific, we have first calculated 
the total cross section a P_ rod > PYTHIA (m H ± ) with the built-in routine for ti production in PYTHIA for 
all m H ± = 80, 100, 130 and 150 GeV and then calculated from this the respective rescaling factors 
c(m H ±) = 5.2 pb/o"^ rod ' PYTHIA (m^±) for each m H ±. Then we have calculated the background cross 
section for m H ± = 80 GeV into the final state with the signature 2b + 2j + Tj et + p™ lss by enforcing 
the respective decay channels in PYTHIA using the built-in routine for tt production and multiplied it 
with c(80 GeV). In the same manner we have calculated the signal cross sections with the PYTHIA 
routines for tbH^ production by enforcing the respective decay channels in PYTHIA and multiplying 
with the rescaling factors c(m H ±) for m H ± = 80, 100, 130, 150 GeV. The resulting cross sections are 
given in Table Q] before (cr th ) and after (a) applying the basic cuts p{ ets > 20 GeV and the hard cut 
^miss y log GeV. For the four signal masses, the tbH^ and tt — > tbH^ cross section calculations agree 
numerically. 



Table 1: Tevatron cross sections of background qq,gg — * tt and signal qq,gg — * tbH for tan/3 = 30 and m H ± = 
80, 100, 130 and 150 GeV into the final state 2b + 2j + r jct + p? iss before (cr th ) and after (a) the basic cuts (p t > 20 GeV for 
all jets) and the hard cut ( P r ss > 100 GeV). 





qq, gg -> tt 


q 


q,gg - 


^tbH± 


m H ± (GeV) 


80 


80 


100 


130 


150 


a th (fb) 


350 


535 


415 


213 


85 


a (fb) for pj ets > 20 GeV 


125 


244 


202 


105 


32 


u (fb) for {ri ets ,pf ss ) > (20, 100) GeV 


21 


30 


25 


18 


7 



2.3 Event Preselection and Discussion of Discriminant Variables 

The expected cross sections of the 2b + 2j + Tj e t + pf lss signature are of the same order of magnitude 
for the signal and background reactions, as shown in Table Q] Thus, the same number of signal and 
background events is assumed for the analysis of different kinematic selection variables. For the signal 
5 • 10 5 events have been simulated with PYTHIA for each charged Higgs mass at the Tevatron energy 
of 1 .96 TeV using the built-in ti routine in the ti — > tbH^ approximation, while for the ti background 
also 5 • 10 5 events have been simulated using the built-in ti routine. Then the basic cuts p} cts > 20 GeV 
are applied. An additional hard cut on the missing transverse momentum p™ lss > 100 GeV is used to 
suppress the QCD background, as for example demonstrated in Ref. ||30l . After the additional anti-QCD 
cut about 28000 to 42000 signal events, depending on the simulated charged Higgs bosons mass, and 
about 30000 it background events remain. Other background reactions, for example W+jet production, 
are expected to be negligible because they have either a much lower production cross section or are 
strongly suppressed compared to background, as quantified for example in Ref. ll30l . In addition to the 
previous study (based on 5000 x BR(t — > hadrons) events each) [33], the present one applies an IDA 
method 071 to explore efficiencies and purities. As already mentioned, particular attention is devoted to 
the study of spin sensitive variables in the exploitation of polarization effects for the separation of signal 
and background events. 

Figures [T|-[8] show examples of the signal and background distributions of some of the kinematic 
variables used in the IDA method and the respective difference between signal and background distribu- 
tions. namely: 

3 Clearly, for a proper experimental study, factorization and renormalization scales for our background process qq, gg — > 
ti — > tbW ± ought to be chosen appropriately, i.e., unrelated to the charged Higgs boson mass. 
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• the transverse momentum of the r jet, p[ set (Fig.Q}, 

• the transverse momentum of the leading 7r ± in the r jet, pf^ (Fig. O 

• the ratio pf jp^ (Fig.©, 

• the transverse momentum of the second (least energetic) b quark jet, p\ 2 (Fig. HJ), 

• the transverse mass in the Tj et + p™ lss system, mt = \j2pt° t pf Llss [l — cos(A</>)], where A(p is the 
azimuthal angle between p T t sct and p™ lss (Fig. [50 

• the invariant mass distribution of the two light quark jets and the second b quark jet, mjjb 2 (Fig. [6]), 

• the spatial distance between the r jet and the second b quark jet, AR(t, 62) = \J (Ac/)) 2 + (At]) 2 , 
where Ac/) is the azimuthal angle between the r and b jet (Fig. [7]) and 

• the sum of the (scalar) transverse momenta of all the quark jets, iJj e t s = pf + pf + p\ 1 + v\ z 
(Fig.©. 

The distributions of signal and background events are normalized to the same number of 10 4 events, in 
order to make small differences better visible. 

The signal and background distributions for the variables shown in Figs. [4]-[8]are as expected rather 
similar for m H ± = m w ± and are hence mostly important to discriminate between signal and background 
in the IDA for m H ± > m w ±. Especially the transverse mass, Fig. [51 shows a large variation with the 
charged Higgs boson mass. However, the different spin of the charged Higgs boson and the W boson 
has a large effect on the r jet variables p T t sct and pj (Figs. Q] and O resulting in significantly different 
distributions of signal and background even for m H ± = m w ± . Moreover, the spin effects in the p t iat and 
pi distributions are correlated which can be seen in Fig. [3] where the distributions of the ratio pi /p T t >et 
|[26l l28l |29l show even larger differences. This highlights the importance of the additional variable p\ 
(and hence pi /pl' ct ), compared to a previous study (33j. The large separation power of this variable is 
indeed due to the different r polarizations in signal and background as can be inferred from the lower 
plots in Figs. [l|-[3j There the signal and background distributions for p T t sct , pi and pi /p T t sct are shown 
for reference samples where the r decay has been performed without the inclusion of spin effects with 
the built-in routines of PYTHIA and hence the differences between signal and background nearly vanish. 

2.4 Iterative Discriminant Analysis (IDA) 

The IDA method is a modified Fisher Discriminant Analysis [ 37 ] and is characterized by the use of a 
quadratic, instead of a linear, discriminant function and also involves iterations in order to enhance the 
separation between signal and background. 

In order to analyze our events with the IDA method, signal and background have been split in two 
samples of equal size. With the first set of samples the IDA training has been performed and then the 
second set of samples has been analyzed. We have used the following 20 variables in the IDA study: 
the transverse momenta p[ >ct , pi , p™ lss , p^ 1 , p\ 2 , p 3 ^ , p 32 , pi 3 ; the transverse mass m t ; the invariant 
masses rrijj, rrijjbj, rrijj^, and s = rrijj^; the spatial distances AR(t, 61), AR(t, 62), AR(r,ji), 
AR(t, J2); the total transverse momenta of all quark jets iJj e ts and of all jets H^w = Hj ets + p]' ct . In 
the analysis of real data, b-quark tagging probabilities and the reconstruction of t and W masses could 
be used to improve the jet pairing, and replace the allocation of least and most energetic 6-jet by a 
probabilistic analysis. 

The results of the IDA study are shown in Figs. [9] and [10] for the event samples with spin effect 
in the r decays for m H ± = 80, 100, 130, 150 GeV and for the reference samples without the spin effect 
for mjj± = 80 GeV in order to illustrate the spin effect. In all plots of the IDA output variable the 
number of background events has been normalized to the number of signal events. Two IDA steps have 
been performed. Figure [9] shows the IDA output variable after the first step, where 90% of the signal is 
retained when a cut at zero is applied. The signal and background events after this cut are then passed 
to the second IDA step. Figure [10] shows the IDA output variable distributions after the second step. 
A cut on these distributions leads to the efficiency and purity (defined as ratio of the number of signal 
events divided by the sum of signal and background events) combinations as shown in the lower right 

4 Strictly speaking this is not the transverse mass since there are two neutrinos in the decay chain of the charged Higgs boson 
we are considering, even so the characteristics of this mass are very similar to that of the true transverse mass. 
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plot in Fig. [TO] These combinations define the working point (number of expected background events 
for a given signal efficiency) and the latter can be optimized to maximize the discovery potential. The 
difference between the dashed (no spin effects in r decay) and solid (with spin effects in r decay) lines 
for rriH± = 80 GeV in the lower right plot in Fig. [10] stresses again the importance of the spin effects to 
separate signal and background. 

In order to illustrate the effect of the hard cut on the missing transverse momentum (pf lss > 
100 GeV), which is imposed to suppress the QCD background, the final efficiency-purity plot of the 
IDA analysis is shown in Fig. (TTJ for m H ± = 80 GeV for two reference samples (red, long dashed: 
with spin effects in the r decay; red, dotted: without spin effects) without imposing the hard cut. The 
black lines (dashed and solid) are for the samples with the hard cut as also shown in the lower right 
plot in Fig. |T0j As expected the achievable purity for a given efficiency decreases with the hard cut, 
therefore the spin effects become even more important to separate signal and background. In principle, 
by choosing the signal reduction rates in the previous IDA iterations, the signal and background rates in 
the final distributions can be varied appropriately. However, we have checked that a different number of 
IDA iterations and/or different efficiencies for the first IDA iteration have only a minor effect on the final 
result. 
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Fig. 1: p t distributions of the r jet for the tbH ± signal and the tt background for yfs = 1.96 TeV (left) and the respective 
differences between signal and background (right). The lower plots show distributions without spin effects in the r decays. 
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Fig. 2: p t distributions of the leading from the r decay for the tbH ± signal and the ti background for yfs — 1.96 TeV 
(left) and the respective differences between signal and background (right). The lower plots show distributions without spin 
effects in the r decays. 
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Fig. 4: p t distributions of the second (least energetic) b quark jet for the tbH signal and the tt background for yfs - 
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Fig. 9: Distributions of the IDA output variable in the first IDA step for the tbH signal (solid, red) and the tt background 
(dashed, blue) for ^/s = 1.96 TeV. 
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Fig. 10: Upper row, middle row and lower left figure: distributions of the IDA output variable in the second IDA step for 90% 
efficiency in the first IDA step (corresponding to a cut at in Fig. [9} for the tbH ± signal (solid, red) and the ti background 
(dashed, blue). Lower right figure: efficiency as a function of the purity when not taking the spin effects in the r decay into 
account for m H ± = 80 GeV (dashed) and with spin effects in the r decay for m H ± — 80, 100, 130, 150 GeV (solid, from left 
to right). Results are for the Tevatron. 
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Purity 

Fig. 11: Efficiency as a function of purity for m H ± = 80 GeV and yfs = 1.96 TeV. The black lines are the results after 
applying the hard cut p™" 58 > 100 GeV when not taking the spin effects in the r decay into account (dashed) and with spin 
effects in the r decay (solid), as also shown in FigQJj] The red lines are the results without applying the hard cut on p" llss when 
not taking the spin effects in the r decay into account (dotted) and with spin effects in the r decay (long dashed). 

3. LHC ENERGY 

The simulation procedure and the emulation of the detector response are the same as those outlined in 
Sect. 2.1 for the Tevatron, as well as, for the preselection and IDA method, as described in Sects. 2.3 
and 2.4, respectively. Hence, only the expected LHC rates are discussed, followed by the description of 
changes in the distributions of the variables and the final IDA results. 

Unlike the case of the Tevatron, where only charged Higgs masses smaller than the top quark 
mass can be explored, and 2HDM/MSSM signatures practically rely on tv t pairs only, at the LHC 
the phenomenology is more varied. Here, the search strategies depend strongly on the charged Higgs 
boson mass. If m#± < mt — nib (later referred to as a light Higgs boson), the charged Higgs boson 
can be produced in top (anti-)quark decay. The main source of top (anti-)quarks at the LHC is again 
ti pair production (cr t j = 850 pb at NLO) IT381 . For the whole (tan/3, m H ±) parameter space there is 
a competition between the bW^ and bH^ channels in top decay keeping the sum BR(i — ► bW + ) + 
BR(t — > bH + ) at almost unity. The top quark decay to bW^ is however the dominant mode for most 
of the parameter space. Thus, the best way to search for a (light) charged Higgs boson is by requiring 
that the top quark produced in the tbH^ process decays to a W . While in the case of decays r's 
will be tagged via their hadronic decay producing low-multiplicity narrow jets in the detector, there are 
two different W decays that can be explored. The leptonic signature bbH^W 11 — ► bbrvlv provides 
a clean selection of the signal via the identification of the lepton I = e, \i. In this case the charged 
Higgs transverse mass cannot be reconstructed because of the presence of two neutrinos with different 
origin. In this channel charged Higgs discovery will be determined by the observation of an excess of 
such events over SM expectations through a simple counting experiment. In the case of hadronic decays 
bbH^W^ — > bbrvjj the transverse mass can instead be reconstructed since all neutrinos are arising 
from the charged Higgs boson decay. This allows for an efficient separation of the signal and the main 
ti — ► bbW^W^ — > bbrujj background (assuming m H ± > m w ±). The absence of a lepton (e or /x) 
provides a less clean environment but the use of the transverse mass makes it possible to reach the same 
mass discovery region as in the previous case and also to extract the charged Higgs boson mass. Both 
these channels show that after an integrated luminosity of 30 fb _1 the discovery could be possible up to 
a mass of 150 GeV for all tan/3 values in both ATLAS and CMS Il9ll39ll40l. 

If the charged Higgs is heavier than the top quark, the dominant decay channels are — ► TV 
and -> tb depending on tan j3. They have both been studied by ATLAS and CMS I4fll42"ll43ll44l . 
The charged Higgs bosons are produced in the pp — > tbH^ channel. For the if ± — ► tb decay, a charged 
Higgs boson can be discovered up to high masses (m H ± ~ 400 GeV) in the case of very large tan j3 
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values and this reach cannot be much improved because of the large multi-jet environment. For the 
— > rv decay mode this reach is larger due to a cleaner signal despite a lower BR. In this case the 5a 
reach ranges from tan = 20 for m H ± = 200 GeV to tan j3 = 30 for m H ± = 400 GeV. 

For the LHC, signal and background events have been simulated in the same way as for the Teva- 
tron as described before, however, without implying any rescaling factor to match a measured tt cross 
section. Table [2] lists the resulting cross sections before (a th ) and after (a) applying the basic cuts 
p\ ets > 20 GeV and the hard cut p™ lss > 100 GeV. The LHC rates allow for the discovery to be less 
challenging than at the Tevatron in the region m H ± ~ m w ± , yet the separation of signal events from 
background remains crucial for the measurement of the charged Higgs mass. 

Table 2: LHC cross sections of background qq, gg — > tt and signal qq, gg — > tbH ± for tan /3 = 30 and m H ± = 80, 100, 130 
and 150 GeV into the final state 2b + 2j + r jct + p™ iss before (cr th ) and after (a) the basic cuts (p t > 20 GeV for all jets) and 
the hard cut (p? iaa > 100 GeV). 





qq, gg -»■ tt 


qq, 99 - 


+ tbH± 


m H ± (GeV) 


80 


80 


100 


130 


150 


(7 th (pb) 


45.5 


72.6 


52.0 


24.5 


9.8 


cr (pb) for pf ts > 20 GeV 


17.3 


33.9 


25.7 


12.2 


3.8 


<r (pb) for (p\ ets ,pf ss ) > (20, 100) GeV 


4.6 


6.0 


4.8 


2.9 


1.2 



The kinematic distributions are shown in Figs. [T2]to[T9]for y/s = 14 TeV The choice of variables 
is identical to the one for the Tevatron and allows for a one-to-one comparison, the differences being 
due to a change in CM energy (and, to a somewhat lesser extent, due to the leading partonic mode of 
the production proces^]). The main differences with respect to Figs. [T|-[8] are that the various transverse 
momenta and invariant masses have longer high energy tails. In particular, it should be noted that the 
effect of the spin differences between and events can be explored very effectively also at LHC 
energies, e.g. the ratio pf /p t lct is shown in Fig. [14] which is very sensitive to the spin effects. These 
observations lead to the conclusion that the same method using spin differences can be used to separate 
signal from background at both the Tevatron and the LHC. 

The distributions of the IDA output variables for the study at = 14 TeV for two steps with 
90% efficiency in the first step are shown in Figs. [20] and [21] These distributions are qualitatively similar 
to those for the Tevatron (Figs. 191 and[T0l) and the final achievable purity for a given efficiency is shown in 
Fig. ED As for the Tevatron energy a good separation of signal and background events can be achieved 
with the spin sensitive variables and the IDA method even in case m H ± ~ m w ±. For heavier 
masses the separation of signal and background events increases due to the kinematic differences of the 
event topology. 



5 As the latter is dominated by qq annihilation at the Tevatron and gg fusion at the LHC. 
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Fig. 12: p t distributions of the r jet for the tbH ± signal and the tt background for yfs = 14 TeV (left) and the respective 
differences between signal and background (right). The lower plots show distributions without spin effects in the r decays. 
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Fig. 13: pt distributions of the leading tt* from the r decay for the tbH ± signal and the tt background for y/s = 14 TeV (left) 
and the respective differences between signal and background (right). The lower plots show distributions without spin effects 
in the r decays. 
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Fig. 14: Distributions of the ratio pi /Pt JCt for the tbH ± signal and the ti background for ^fs — 14 TeV (left) and the 
respective differences between signal and background (right). The lower plots show distributions without spin effects in the r 
decays. 
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Fig. 17: Invariant mass distributions of the two light quark jets and the second (least energetic) b quark jet for the tbH signal 
and the tt background for y/s = 14 TeV (left) and the respective differences between signal and background (right). 
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Fig. 18: Spatial distance AR(t, 62) = \J (A<^) 2 + (A77) 2 distributions (where A</> is the azimuthal angle in rad between the 
r and b jet) for the tbH ± signal and the ti background for -y/s = 14 TeV (left) and the respective differences between signal 
and background (right). 
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Fig. 19: Distributions of the total transverse momentum of all quark jets, Hj cts = p^ 1 + p J t 2 + p t 1 + p t 2 , for the tbH signal 
and the ti background for y/s = 14 TeV (left) and the respective differences between signal and background (right). 
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Fig. 20: Distributions of the IDA output variable in the first IDA step for the tbH ± signal (solid, red) and the tt background 
(dashed, blue) for y/s = 14 TeV. 
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Fig. 21: Upper row, middle row and lower left figure: distributions of the IDA output variable in the second IDA step for 90% 
efficiency in the first IDA step (corresponding to a cut at in Fig. 120b for the tbH ± signal (solid, red) and the ti background 
(dashed, blue). Lower right figure: efficiency as a function of the purity when not taking the spin effects in the r decay into 
account for m H ± = 80 GeV (dashed) and with spin effects in the r decay for m H ± — 80, 100, 130, 150 GeV (solid, from left 
to right). Results are for the LHC. 
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4. CONCLUSIONS 



The discovery of charged Higgs bosons would be a clear sign of physics beyond the SM. In this case 
study we have investigated charged Higgs boson topologies produced at the current Tevatron and LHC 
energies and compared them against the irreducible SM background due to top-antitop production and 
decay. While sizable differences between signal and background are expected whenever m#± 7^ m^± , 
near the current mass limit of about m H ± 80 GeV the kinematic spectra are very similar between 
SM decays and those involving charged Higgs bosons. In this case, spin information will significantly 
distinguish between signal and irreducible SM background. In fact, we have considered hadronic tv t 
decays of charged Higgs bosons, wherein the r polarization induced by a decaying (pseudo)scalar object 
is significantly different from those emerging in the vector (W^) decays onsetting in the top-antitop 
case. For a realistic analysis which is not specific for a particular detector, a dedicated Monte Carlo 
event generation and a simplified multipurpose detector response approximation have been applied. The 
identification of a hadronic tau-lepton will be an experimental challenge in an environment with typically 
four jets being present. We have demonstrated how an IDA method can be an applied to separate signal 
and background when the differences between the signal and background distributions are small. Our 
results show that the IDA method will be equally effective at both the Tevatron and LHC. While only 
the dominant irreducible ti background has been dealt with in detail, we have also specifically addressed 
the QCD background. A suitably hard missing transverse momentum cut has been applied to reject 
such jet activity and we have demonstrated that although the discriminative power is reduced by such 
a cut, the reduction is small compared to the gain from including the r polarization effects. Using the 
differences in r polarization between the signal and the dominant SM irreducible ti background is crucial 
for disentangling the former from the latter. 
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